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SIGNALS  IN  NEUROMEDIATE  SYNCOPE 

F.  Censi1,  G.  Calcagnini1,  S.Strano2,  P.  Bartolini1,  V.  Barbara1 

'Biomedical  Engineering  Laboratory,  Istituto  Superiore  di  Sanita’,  Roma,  Italy 
2Clinical  Pathophysiology,  La  Sapienza  University  of  Rome,  Italy 

Abstract-  Aim  of  this  study  is  to  evaluate  the  degree  of  variables  has  indeed  improved  the  knowledge  about  the 

coupling  among  the  cardiovascular  variability  series  and  autonomic  functioning,  even  if  applications  of  such  analyses 

the  respiration  in  subjects  susceptible  to  neurally  to  investigate  autonomic  nervous  system-linked  pathologies 

mediated  syncope,  in  comparison  to  normal  subjects.  21  are  scarce.  In  addition,  given  the  complexity  of  the 

patients  susceptible  to  syncope  and  10  sex  and  age  mechanisms  involved  in  the  autonomic  regulation,  the 

matched  control  subjects  have  been  enrolled  in  the  study.  analysis  should  take  into  account  the  non-linear  features  of 

ECG,  respiration  activity,  and  arterial  blood  pressure  the  coupling  between  the  involved  variables, 

were  simultaneously  recorded  at  rest  (controlled  and  free  Aim  of  this  study  is  to  evaluate  the  degree  of  coupling 

breathing)  and  during  70°  head-up  TILT  test  (free  among  the  cardiovascular  variability  series  and  the 

breathing).  The  degree  of  non-linear  coupling  among  respiration  in  subjects  susceptible  to  neurally  mediated 

heart  rate  variability  (HRV),  blood  pressure  variability  syncope,  in  comparison  to  normal  subjects.  The  degree  of 

(BPY)  and  respiration  has  been  quantified  by  means  of  a  non-linear  coupling  has  been  quantified  by  means  of  a 

multivariate  embedding -based  approach.  11  patients  multivariate  embedding-based  approach,  recently  proposed 

developed  syncope  during  the  TILT  test.  The  non-linear  by  Hoyer  et  al.  [7].  Particularly,  we  used  two  non-linear 

coupling  among  the  cardiovascular  variables  and  coupling  indexes  called  independence  of  complexity  (IC)  and 

respiration  turned  out  to  be  stronger  during  controlled  independence  of  predictability  (IP), 
breathing  than  during  free  breathing,  both  for  TILT- 

positive  and  TILT-negative  group.  During  late  TILT  D.  METHODOLOGY 

phase,  in  the  TILT-positive  group,  the  non-linear 

coupling  showed  a  significant  increase.  If  the  proposed  A.  Experimental  protocol  and  study  population 
non-linear  coupling  indexes  can  be  considered  expression 

of  the  coupling  mechanisms  involved  in  the  vagal  ECG,  respiration  activity  (RA),  and  arterial  blood 

regulation  of  the  cardiovascular  system,  an  increase  of  the  pressure  (ABP)  were  simultaneously  recorded  at  rest  and 

vagal  tone  accompanied  by  a  decrease  of  the  sympathetic  during  70°  head-up  tilt  test  (figure  1).  Surface  ECG  (II  lead) 

activity  seem  to  occur  before  a  vasovagal  event.  was  recorded  by  an  analog  electrocardiograph  (MCR  I, 

Keywords  -  Non-linear  coupling,  heart  rate  variability,  Esaote,  Italy).  APB  was  continuously  and  non-invasively 

cardiorespiratory  coordination,  neurally  mediated  syncope  recorded  by  photopletismographic  technique  (Finapres, 

Omheda,  USA),  and  RA  was  monitored  using  a 

I.  INTRODUCTION  pletismographic  thoracic  belt.  The  signals  were  sampled  in 

real-time  (sampling  frequency:  500  Hz,  resolution:  12  bit, 
Cardiovascular  variables  are  characterized  by  DT2801A,  Data  Translation,  USA)  and  stored  in  a  magneto 
pseudoperiodic  oscillations  reflecting  the  complex  control  optical  disk  for  further  analysis.  Data  were  collected  in  the 

mechanisms  mediated  by  the  autonomic  nervous  system.  Clinical  Pathophysiology  Unit  of  the  University  of  Rome,  ‘La 

Over  the  last  20  years  a  great  deal  of  research  regarding  Sapienza’.  All  experiment  were  done  in  a  light  attenuated 

short-term  heart  rate  and  blood  pressure  variability  has  been  room  starting  from  2  PM 

carried  out  in  frequency  domain  [1][2][3].  This  approach 
significantly  helped  to  better  understand  the  physiological 
meaning  of  the  two  main  rhythms  observed  in  heart  rate  and 
blood  pressure  fluctuations:  a  low  frequency  component  (LF, 

0.04-0.15  Hz)  which  is  believed  to  be  an  expression  of  the 
baroreflex  control,  mediated  by  both  sympathetic  and 
parasympathetic  branches,  and  a  high  frequency  component 
(HF,  0.15-0.4  Hz)  synchronous  with  respiration,  generally 
considered  a  marker  of  the  parasympathetic  activity  [4]. 

In  parallel  with  the  above  linear  approaches,  the 

investigation  of  non-linear  dynamics  in  heart  rate  variability 
have  improved  our  understandings  regarding  the 

physiological  and  pathologically  disturbed  behaviors  of  the 
autonomic  nervous  system  [5][6].  Fig.  1.  Experimental  setup. 

It  should  be  kept  in  mind,  however,  that  other 

physiological  variables  as  well  as  external  stimuli  are  In  supine  position,  two  5-minute  stages  of  controlled 

involved  in  the  autonomic  nervous  system-mediated  respiration  have  been  performed  (8  and  12  breaths/minute, 

regulation.  The  analysis  of  the  interaction  among  these  0.13  Hz  and  0.20  Hz  respectively),  followed  by  a  5-minute 
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stage  of  spontaneous  respiration  (REST,  table  I).  Upright 
posture  was  maintained  until  the  development  of  syncope  or 
for  a  maximum  duration  of  40  minutes.  Head-up  tilt  was 
defined  positive  if  hypotension  (systolic  blood  pressure  lower 
than  80  mmHg),  accompanied  by  slowing  of  the  heart  rate 
(lower  than  40  bpm)  developed.  Subjects  with  syncope  were 
quickly  returned  to  the  horizontal  position.  After  return  of 
patients  to  supine  posture,  a  5minute  recovery  stage  has  been 
considered. 

TABLE  I 

STAGES  OF  THE  EXPERIMENTAL  PROTOCOL 


Stage 

breaths/minute 

Duration 

(minutes) 

1 

REST  -  supine  position 

12 

5 

Instrumentation 

free 

3 

calibration 

2 

8 

5 

Instrumentation 

free 

3 

calibration 

3 

free 

5 

TILT  -  70°  head  up  tilt  test 

4 

free 

<40 

RECOVERY  -  supine  position 

5 

free 

5 

The  study  population  consisted  of  21  patients  (10  women) 
with  two  or  more  episodes  of  unexplained  syncope  during  the 
preceding  6  months,  and  of  10  sex  and  age  matched  control 
subjects.  No  patients  had  evidence  of  structural  heart  disease 
or  diabetes.  Patients'  average  age  was  24±7  years  (range: 
14444). 

B.  Data  pre-processing 

The  low-pass  filtered  event  series  (ES)  was  used  to 
extract  the  heart  rate  variability  (HRV)  signal  since  it  is  a 
reliable  time  domain  representation  with  high  temporal 
resolution.  In  ES  each  beat  is  replaced  by  a  8  function;  the 
signal  can  be  described  as: 

N- 1 

x(t)=  £S(f-f  k) 

k= 0 

where  8  is  the  Dirac  delta  function,  f  is  the  occurrence  time 
of  the  It1  beat  and  N  is  the  total  number  of  beats.  To  obtain 
the  heart  rate  variability  series,  the  ES  was  low-pass  filtered 
at  0.5  Hz  and  re-sampled  at  2  Hz  [8].  Correspondingly,  blood 
pressure  variability  (BPV)  and  RA  series  were  obtained  by  a 
similar  low-pass  filtering  procedure  (FIR,  10  order,  cut-off 
frequency  0.5  Hz)  and  a  re-sampling  at  2  Hz  of  the  ABP  and 
RA,  respectively.  According  to  this  procedure,  HRV  and  RA 
series  were  expressed  in  arbitrary  units  (a.u.),  and  BPV  in 
mmHg.  This  approach  particularly  guarantees  the  time 
synchronization  between  HRV,  BPV  and  RA. 

C.  Non-linear  coupling  quantification 

The  non-linear  coupling  quantification  is  based  on  a 
multivariate  embedding  approach.  The  multivariate 
embedding  procedure  allows  for  the  reconstruction  of  the 
joint  behaviour  of  multiple  systems  which  are  more  or  less 
dependent  [5][9].  According  to  the  Takens’  reconstruction 
theorem  [10],  the  dynamics  of  a  system  can  be  reconstructed 
by  taking  the  time-delayed  values  of  one  single  observed  time 


series.  The  joint  system  dynamics  of  two  coupled  subsystems 
(x(t)  and  y(t))  can  be  reconstructed  by  an  analogous  delay 
vector  -  q(t)  -  using  values  of  the  two  contributing 
subsystems,  as  follows: 

q(t)=  [x(t),  x(t-x),  x(t-2x),  ...  ,  x(t-nrx),  y(t),  y(t-x),  y(t-2x), 
•  ••,  y(t-nx)] 

where  x(t)  and  y(t)  have  embedding  dimension  m  and  n 
respectively,  and  the  joint  embedded  system  has  dimension 
m+n. 

The  method  consist  on  the  computation  of  the  correlation 
dimension  (CD)  and  of  the  correlation  entropy  (CE)  for  the 
individually  embedded  series  (x(t)  and  y(t))  as  well  as  for  the 
jointly  embedded  series  (q(t)),  and  on  the  comparison  of  the 
obtained  values.  CD  and  CE  have  been  estimated  according 
to  the  Grassberger-Procaccia  method  for  the  estimation  of  the 
correlation  integral  [11].  The  correlation  integral  counts  for 
the  number  of  pairs  of  points  x(ti),x(tj)  which  are  separated  by 
a  distance  less  than  e  in  a  phase-space  of  embedding 
dimension  m,  according  to  the  following  formula: 

C  (e,  m)  =  — J^H(e-\x(t.)-  x(  t. )  I) 

N  ;= i  j=i+ 1 

where  N’=N-m+l,  N  is  the  number  of  points  of  the  time 
series,  and  H  is  the  Heaviside  step  function  (H(y)=0  for  \<0, 
and  H(y)=l  for  y>0). 

We  estimated  CD  as  the  slope  of  the  log-log  plot  of 
C(£,m)  vs.  ein  the  scaling  region,  as  follows  [9]: 


CD  (e) 

m v  7 


d  log  C  (e,  rri) 
fifloge 


In  the  scaling  region  CD  is  independent  of  e  and  m.  Since 
in  physiological  data  sets  no  convincing  scaling  region  exists, 
we  performed  the  comparison  of  different  time  series  by 
using  appropriate  standard  values  of  e  and  m  [5].  CE  was 
estimated  as: 


CEJt)  =  -  In 
x 


C(e,m) 
C(£,m  + 1) 


where  eand  m  were  chosen  according  to  the  CD  estimation. 

IC  has  been  extracted  by  combining  the  CDs  of  each 
individually  embedded  series  (CDX  and  CDy)  with  that  of  the 
jointly  embedded  signal  (CDq),  as  follows  [7]: 


I  CD  -CD  I  + 1  CD  —CD  I 
]r  —  i _ * _ l _ >_ 

CDX  +  CD,, 

Similarly,  IP  has  been  extracted  by  combining  the  CE 
values  as  follows  [7]: 


,  ,,  I  CEq  -  CEX I  + 1  CEq  -  CEy  I 
CEx  +  CEy 

Both  indexes  vary  between  0  and  1.  0  indicates  fully 
coupled  quantities  and  1  stands  for  completely  independent 
variables.  By  means  of  these  indexes,  we  estimated  the  non¬ 
linear  coupling  between  HRV  and  respiration,  BPV  and 
respiration,  and  HRV  and  BPV. 

In  this  work,  approximate  (coarse-grain)  correlation 
dimension  and  correlation  entropy  have  been  calculated  for 
an  embedding  dimension  ;«=10  (m= 20  for  the  joint  embedded 
vector)  and  a  distance  r  equal  to  10%  the  maximum  distance 
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of  points  in  the  phase  space.  The  time  delay  x  has  been 
chosen  according  to  the  first  zero  of  the  auto-correlation 
function  of  each  series,  which  for  the  analyzed  series 
corresponds  to  about  one  quarter  of  the  characteristic  period 
(x=10). 

m.  RESULTS 

11  out  of  21  subjects  developed  syncope  during  head-up 
TILT  test.  The  TILT  -positive  subjects  developed  pallor, 
diaphoresis,  dizziness  and/or  transient  loss  of  consciousness 
within  5  to  15  minutes  of  head-up  TILT.  The  TILT  was 
terminated  when  systolic  blood  pressure  fell  below  80  mmHg 
or  when  HR  fell  below  40  bpm.  All  subjects  recovered 

spontaneously  on  return  to  the  supine  position. 

The  coupling  between  respiration,  HRV  and  BPV  has 
been  analyzed  during  the  two  controlled  breathing  stages  (5 
minutes  each),  during  REST  (5  minutes),  TILT,  and  the 
recovery  phase.  Because  the  period  of  upright  TILT  varied  in 
duration  among  patients  with  syncope  depending  on  whether 
symptoms  developed  and  the  time  that  syncope  occurred, 

total  duration  of  upright  posture  in  each  patient  or  subject  was 
divided  into  3  intervals  indicated  as:  early  TILT,  middle 
TILT  and  late  TILT  (which  corresponds  to  pre-syncope  in 

TILT  positive  group).  The  recovery  phase  started  after  return 

of  the  patient  to  supine  posture  and  lasted  5  minutes  (figure 
2).  Linear  coupling  has  been  quantified  as  the  maximum  of 
the  cross-correlation  function. 

Figure  3  shows  the  values  of  the  coupling  indexes  (IC,  IP 
and  cross-correlation)  between  HRV  and  respiration 
(HRV/respiration)  in  the  7  phases  of  the  protocol,  for  the 
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Fig.  2.  Phases  of  the  experimental  protocol. 


TILT -positive  subjects  (left  panel)  and  for  the  TILT-negative 
group  (right  panel)  expressed  as  mean  ±  standard  deviation. 
The  results  obtained  for  the  TILT -negative  group  are  not 
significantly  different  from  those  obtained  for  the  control 
group. 

The  non-linear  coupling  indexes  show  that  the  coupling  is 
significandy  stronger  during  the  8  breaths/minute  stage  than 
during  12  breathsAmnute  stage  for  both  TILT -positive  and 
TILT-negative  group  (p<0.01).  In  addition,  during  controlled 
respiration,  the  coupling  is  significantly  higher  than  during 
free  breathing  at  REST  for  both  groups  (p<0.01).  From  REST 
condition  to  early  and  middle  TILT,  the  coupling  does  exhibit 
a  significant  decrease  for  both  groups  for  HRV/respiration, 
BPV/respiration  (p<0.01).  During  late  TILT  phase,  in  the 
TILT -positive  group  the  non-linear  coupling  has  a  significant 
increase  (p<0.01)  and  during  the  recovery  phase  i  returns  to 
values  similar  to  that  observed  during  the  REST  phase. 
Linear  coupling  does  not  significantly  vary  in  the  TILT 
phases. 

Similar  results  have  been  obtained  for  the  coupling 


HRV/respiration 


IC 


IP 


Linear 

coupling 


TILT  -positive  subjects 


Fig.  3 .  Coupling  between  HRV  and  respiration:  non-linear  indexes  values  during  the  phases  of  the  controlled  breathing  and  TILT  -test  protocol,  for  the 
TILT  -positive  subjects  (left  panel)  and  the  TILT -negative  group  (right  panel).  Dashed  line  denotes  the  increase  of  the  coupling  strenght  during  the  late 

TILT  phase  (or  pre-syncope)  for  the  TILT -positive  subjects. 
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BPV/respiration  and  HRV/BPV  (not  shown). 

No  significant  difference  has  been  observed  between  TTLT- 
positive  and  TILT-negative  group  during  controlled 
breathing,  and  free  breathing  at  REST  and  during  early  and 
middle  TILT.  Of  course  the  difference  does  reach  a  statistical 
significance  during  late  TILT  (p<0.01). 

The  LF/HF  ratio,  computed  during  the  free  breathing  REST 
phase  and  during  the  3  TILT  phases,  was  not  different 
between  the  2  groups  in  all  phases.  During  REST  its  value  is 
significantly  different  from  that  obtained  in  the  early  TILT, 
for  both  groups  (p<0.01).  However,  during  early  TILT  phase 
the  values  obtained  for  LF/HF  ratio  widely  vary  among 
patients.  During  late  TILT,  the  LF/HF  ratio  is  not  different 
from  the  preceding  TILT  phases. 

IV.  DISCUSSION 

The  results  obtained  from  this  investigation  show  an 
increase  of  the  non-linear  coupling  between  cardiovascular 
variability  signal  and  respiration  in  concomitance  with  a 
syncope  event,  which  cannot  be  detected  by  traditional  linear 
methods.  Since  the  autonomic  nervous  system,  governing  the 
links  between  heart  rate,  blood  pressure  and  respiration,  is  a 
complex  system,  it  is  not  surprising  that  a  deeper  non-linear 
analysis  of  the  coupling  among  the  contributing  subsystems 
provides  significant  results. 

Since  respiration  is  usually  closely  coupled  to  heart  rate, 
the  cardiorespiratory  coupling  could  be  of  paramount 
importance  for  the  understanding  of  the  underlying 
pathogenic  mechanisms,  as  well  as  for  the  choice  of  the 
proper  therapeutic  approach.  During  controlled  and  free 
breathing  at  REST,  the  cardiorespiratory  coupling 
(HRV/respiration  and  BPV/respiration)  did  not  differ 
between  TILT -positive  and  TILT-negative  subjects.  This 
finding  is  consistent  with  the  previous  results  obtained  by 
Lipsitz  et  al.  [12],  and  indicates  that  the  baseline 
cardiorespiratory  coordination  is  not  altered  in  subjects 
susceptible  to  syncope.  Similar  results  have  been  found 
during  early  and  middle  TILT,  suggesting  an  adaptation  of 
the  cardiorespiratory  dynamics  at  the  beginning  of  the  head- 
up  TILT  test.  In  addition,  during  the  first  2  TILT  phases,  the 
non-linear  coupling  between  both  HRV  and  respiration  and 
BPV  and  respiration  resulted  to  be  significantly  lower  than 
that  observed  during  REST.  This  finding  is  somehow 

consistent  with  the  de-coupling  found  by  Lipsitz  et  al.  3 

minutes  before  syncope,  and  suggests  a  dissociation  between 
respiration  and  cardiovagal  activity  during  the  first  2  TILT 
phases  probably  associated  with  the  increase  of  the 
sympathetic  tone  [2]. 

One  possible  explanation  for  the  dissociation  between 
respiration  and  heart  rate  variability  could  be  the  high 

sympathetic  tone  during  tilt  that  suppresses  the  respiratory 
sinus  arrhythmia  until  sympathetic  withdrawal  before 
syncope  allows  respiratory  sinus  arrhythmia  to  increase. 
Since  the  non-linear  coupling  between  HRV  and  respiration 
is  likely  considered  the  expression  of  the  parasympathetic 
outflow  to  the  heart,  the  progressive  decrease  of  the  non¬ 

linear  coupling  from  REST  to  TILT  could  be  the  expression 


of  the  shift  in  the  balance  between  sympathetic  and 

parasympathetic  activity.  Particularly,  if  the  proposed  non¬ 
linear  coupling  indexes  can  be  considered  expression  of  the 
coupling  mechanisms  involved  in  the  vagal  regulation  of  the 
cardiovascular  system,  an  increase  of  the  vagal  tone 

accompanied  by  a  decrease  of  the  sympathetic  activity  seem 
to  occur  before  a  vasovagal  event. 
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